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Abstract. We study Schrodinger operators on the real line whose potentials 
are generated by an underlying ergodic subshift over a finite alphabet and a 
rule that replaces symbols by compactly supported potential pieces. We first 
develop the standard theory that shows that the spectrum and the spectral 
type are almost surely constant, and that identifies the almost sure absolutely 
continuous spectrum with the essential closure of the set of energies with van- 
ishing Lyapunov exponent. Using results of Damanik-Lenz and Klassert-Lenz- 
Stollmann, we also show that the spectrum is a Cantor set of zero Lebesgue 
measure if the subhift satisfies the Boshernitzan condition and the potentials 
are aperiodic and irreducible. We then study the case of the Fibonacci subshift 
in detail and prove results for the local Hausdorff dimension of the spectrum 
at a given energy in terms of the value of the associated Fricke-Vogt invari- 
ant. These results are elucidated for some simple choices of the local potential 
pieces, such as piecewise constant ones and local point interactions. In the lat- 
ter special case, our results explain the occurrence of so-called pseudo bands, 
which have been pointed out in the physics literature. 



1. Introduction 

Operators associated with strictly ergodic subhifts over finite alphabets have 
been studied in numerous papers since the 1980's. Many of the developments were 
surveyed in [13] . Until quite recently, most of the effort was devoted to discrete 
Schrodinger operators arising in this way. On the other hand, there have been 
several works considering continuum Schrodinger operators [32] , Jacobi matrices 
, and CMV matrices [22 [23J OH] . Along the way it has been realized that many 
tools and results are quite universal in these various instances. On the other hand, 
when implementing these tools and proving these results, one does have to deal with 
some model-dependent features. In other words, one may in general be guided by 
analogies but there does not seem to be an automatic way of carrying over results 
from one model to another. 

The central example in the discrete Schrodinger setting is given by the Fibonacci 
Hamiltonian. Indeed, the foundational papers by Kohmoto, Kadanoff, Tang [33] 
and Siito [39 , 40J studied this particular model and proved for it the spectral prop- 
erties, zero-measure spectrum and singular continuous spectral measures, which 
turned out to be characteristic for the entire class of models. Moreover, the study 
of this key example is additionally motivated by its relevance to the study of qua- 
sicrystals. We refer the reader to the recent survey [14] . which further explores this 
connection. Continuum analogues of the Fibonacci Hamiltonian, especially those 
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with local point interactions leading to the so-called Kronig-Penney model, have 
been studied in [31 HH [23 E01 EH HH S3 - 

The general theory has not yet been worked out to the natural extent possible 
in the continuum Schrodinger setting. Specifically, the paper |32j does not address 
zero-measure Cantor spectrum, and in fact it replaces the use of Kotani theory with 
Remling's non-ergodic version of it. While this generates a more general result on 
singular spectrum, it fails to provide the starting point for a proof of zero-measure 
Cantor spectrum. 

One main motivation for writing this paper is to study the continuum setting 
with the help of Kotani theory in order to prove zero-measure Cantor spectrum in 
the same generality as it is known to hold in for discrete Schrodinger operators. To 
the best of our knowledge, this paper gives the first explicit example of a continuum 
Schrodinger operator whose spectrum is a Cantor set of zero Lebesgue measure^ 
In the discrete case, zero-measure spectrum is known to hold also for the critical 
almost Mathieu operator [21 [35]. This, however, is a very unstable phenomenon 
and quite atypical in the class of operators the almost Mathieu operator belongs 
to. In particular, it is doubtful whether there is a smooth quasi-periodic continuum 
potential so that the associated Schrodinger operator in L 2 (R) has zero-measure 
spectrum. In the discrete setting, zero-measure Cantor spectrum is also typical for 
limit-periodic potentials in the sense of topological genericity; see pQ. No contin- 
uum analog is known, but may possibly hold. The zero-measure Cantor spectrum 
phenomenon is unstable in this scenario as well due to the fact that the periodic 
potentials are dense in the limit-periodic ones. On the other hand, the zero-measure 
Cantor spectrum of discrete Schrodinger operators associated with suitable strictly 
ergodic subshifts is both typical and stable within this class of operators. Thus, the 
latter result is the primary candidate for an attempt to carry over the phenomenon 
of zero-measure Cantor spectrum to the continuum, and this is precisely what we 
accomplish here. 

The second main motivation for writing this paper is to study the spectrum of 
the continuum version of the Fibonacci Hamiltonian, which is of interest for several 
reasons. Since the continuum operator is unbounded, its spectrum is unbounded 
as well, and hence there are potentially interesting questions one can ask about the 
high-energy regime. Related to this, the energy-dependence of the local structure of 
the spectrum may be reduced to the value the so-called Fricke-Vogt invariant takes 
at the energy in question. This value happens to be constant in the discrete case, 
while it is not constant (and hence carries more information) in the continuum 
setting. Secondly, we are able to address a phenomenon that has been pointed 
out in the physics literature. Namely, the Fibonacci Kronig-Penney model may 
exhibit so-called pseudo bands, where the spectrum is locally quite thick and in 
numerical studies (see [31131]) it is a priori not quite clear if the Cantor structure of 
the spectrum breaks down at such energies. We show that this is due to the local 
Hausdorff dimension being one at such points, which in turn comes from zeros of 
the Fricke-Vogt invariant inside the spectrum. The latter phenomenon is impossible 
in the discrete case! 



*It appears to be even the first example of a continuum Schrodinger operator with potential 
in ij 2 oc jf (R) whose spectrum has zero Lebesgue measure. 
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The organization of the paper is as follows. In Section [2] we describe the general 
setting we will work in. That is, we recall the definition of a subshift over a finite 
alphabet and describe how we associate continuum Schrodinger operators. We also 
prove the standard results that, given an ergodic measure on the subshift, the 
spectrum and the spectral type do almost surely not depend on the element of 
the subshift, and the almost sure spectrum does not contain any isolated points. 
In Section [3j we introduce the Lyapunov exponent. To be precise, we introduce 
both Lyapunov exponents that are natural in this setting - one associated with 
the subhift dynamics and one associated with the space variable on the real line - 
and prove that they are multiples of one another. Then we prove another standard 
result, namely that the almost sure absolutely continuous spectrum is given by the 
essential closure of the set of energies where the Lyapunov exponent vanishes. The 
work of Klassert, Lenz, and Stollmann [35] is then used in Section [4] to derive the 
absence of absolutely continuous spectrum whenever the potentials are aperiodic 
and irreducible. Our main point here is that whenever this result applies, the 
Lyapunov exponent is almost everywhere positive. This in turn then implies that 
the spectrum has zero Lebesgue measure, provided that the subshift satisfies the 
Boshernitzan condition; see Section[5j The Cantor structure is a consequence of this 
since the spectrum is always closed and does not contain isolated points as shown 
earlier. This completes the general theory for these continuum operators we wish 
to present. Moving to a special case, we then consider the subshift generated by 
the Fibonacci substitution in Section [U We establish the usual connection between 
the spectrum and the dynamics of the trace map, prove a result expressing the 
local Hausdorff dimension of the spectrum in terms of the value of the Fricke-Vogt 
invariant at the energy in question, and study the latter quantity for two special 
cases, piecewise constant potentials and the Kronig-Penney model. 



2. Subshifts Over Finite Alphabets and Associated Continuum 
Schrodinger Operators 

In this section we introduce the operators we will study in this work and prove 
some general results concerning the spectrum and the spectral type of these oper- 
ators. 

First we introduce the notion of the concatenation to make rigorous our notion 
of "piecing together functions." To that end, assume that for each n € Z, we have 
l n > and a function f n : [0, i n ) — > E. Assume further that J2n>o 4 = E n <o ^« = 
oo (this ensures that the pieced-together function will be defined on all of R) . We 
define the concatenation of the /„ 's by 



fn (x - YJk=l 4J n > and x <E 
fn(x + J2k=n 4 ) n < and x <E 



V P, V n P 

~ Efc=n 4i ~ J2k=n+l 4 



We remark that the empty sums which occur for n = and n — — 1 are defined 
to be equal to by convention. We will notate this concatenation of /„'s by 
/ = ■ ■ |/_ 2 |/_i| /o I/1I/2I ' ' ' ) ■ Note the use of a box to indicate the position 



of the origin in the concatenation. One can also concatenate only finitely many 
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functions - for m < < n, we can produce (^f m \ ■ • • /o ' ■ • \fn) in the same way, 



but it will only be denned on 



( ' ' ' fu-2\fu-l I fujQ l/wi 1/^2 ' ' ' J ' 



Let A be a finite set, called the alphabet. Consider the discrete topology on A 
and the product topology on A z . The shift transformation T : A z — > A z is given 
by (Tlu)(u) = lu(u + 1). Clearly, T is a homeomorphism. A subset ft C _4 Z is called 
T-invariant if T(f2) = f2. Any closed T-invariant set 57 C „4 Z is called a subshift. 

For every a € .A, we choose a length l a > and a real- valued function / a 6 
T 2 (0,£ o ). Let a subshift £1 C A 1 " be given. For cj G Q, we define the Schrodinger 

operator = + V u in T 2 (]R) via V u = 

The potentials Kj so defined belong to Lf ocunii (R) and hence define self-adjoint 
operators in L 2 (M) in a standard way. We will consider the operators {iTjjwesi as 
a family and be interested in statements about H u that hold for many, most, or all 

Let St denote the shift operator (Stf)(x) = f(x — t) on L 2 (R). For ui G fi, let 
U u = Si . One readily sees that U u is unitary and a computation reveals that 
Htu = U^H^U^j. Thus, spectral data will be T-invariant, so one should expect the 
spectrum and spectral parts of to be nonrandom. Indeed, one has the following 
fact: 



Proposition 2.1. Let (£l,T) be a subshift endowed with an ergodic measure \x. 
Suppose A,£ a , f a ,V u , and H u are defined as above. Then there are closed sets 
E, E ac , E sc , E pp C K such that o,{H u ) = X. for • G { , ac, sc, pp} for fi almost 
every lu G f2. 

Proof. The proof is an adaptation of ideas in the literature; see, e.g., [TX)1 fTTl I5T] , 
We will begin by showing that a suitable family of projections is weakly measurable. 

For a fixed interval / C R, let V U {I) = Xi{H u ). Fix f,g G L 2 (R). We shall 
show that the map lu m- {f,Vu(I)g) is measurable. First, for each n G Z + , let 
fiW = ^{-n,...,n} and ^(n) . q _^ Q( n ) the natura l restriction map. For w G f2, 

we consider the cut off potential Vu , given by — I4;Xj(n)> with 1^ — 

(—(i u _ n + • • • + ^u_i)i^wo + ■ ■ • + 4;„)- We see that u n- (/, e ltv ^ n> g) is measurable 
for all n G N, t G K. But then (f,e ltVu 'g} is the limit of measurable functions and 
hence is itself a measurable function of lu. 

By the Trotter product formula, it follows that e ttHul is also a weakly measurable 
funtion of lu. We can approximate characteristic funtions of intervals with trigono- 
metric polynomials to then see that Vu, (I) is a weakly measurable function of lu for 
all intervals /. 

Since Vtu{I) = U*V u {I)Uu, it follows that lu h> tr(7 , UJ (/)) is measurable and 
satisfies tr(7 , Tw(^)) = tr(V u (I)), so ti(V u {I)) is almost surely constant. 

For each pair of rational numbers p < q, let d{p, q) denote the almost sure value 
of ti(V w ((p, q))) and Pi9 the (full measure) set containing those lu G fi for which 
tr("P w ((p, (?))) = d(p,q). Since the set of rational pairs is countable, it follows that 

= C\ P <q ^p,q i s a mu measure subset of f2. 

Suppose E £ a(H Ul ) for some lo% G f2<j- Choose rational p < q so that B G 
(p, g) C p(H Ut ). Then g) = ti(P ul ((p,q))) = 0. For any other lu 2 G f2 , one 
has ti(V UJ2 ((p, <?))) = d(p,q) — 0, so £ ^ o-{H^ 2 ). By symmetry, there is some 
nonrandom closed set ScK with a(H UJ ) = E for all w e O . 
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For the almost sure constancy of spectral parts, let V^f denote orthogonal pro- 
jection onto the absolutely continuous subspace of and define V^f and "PP p anal- 
ogously. Next, let V'(I) — VuiJYPlj, for • e {ac, sc, pp}. Running the argument 
above with V W (I) replaced by V'(I) (and a measurability proof as in [TTT1I3I]) gives 
us a full measure set f2' C f2 and nonrandom sets E ac , E sc , E pp so that cr # (_ff w ) = E, 
for all lu E il' and • G {ac, sc, pp}. □ 

In the usual setting of ergodic dynamically defined potentials, the spectrum lacks 
isolated points. This fact carries over to the present setting. 

Lemma 2.2. Let (O, T) be a sub shift endowed with an ergodic measure [i. Suppose 
•A, £ a , f a ,Vu), H u , and E are defined as above. Then E has no isolated points. 

Proof. It is convenient for us to work with measures on the space of potentials, so 
we let v be the push- forward of /i under the map $ : u M- V u . We also will find 
it convenient to work with IR-invariant families of potentials, so, to that end, let 
f2 = ft x [0, 1). For each (w, t) £ f2, we define V u> t(x) = V u (x + t£ Uo ). Let fx denote 
the product of \x and Lebesgue measure, and v be the pushforward of ft under the 
map * : (u,t) H> V u ,f Denote X = {V u j : ui £ Cl,t £ [0, 1)} 

As in the proof of the previous proposition, we fix an interval I and take V w> t = 
"Pijj,t{I) = Xl{Hu,t) with H u ,t — —A + V Ut t- An argument similar to the one in the 
above proposition establishes that tr(V u ,t) = tr(7 3 w ). 

Let {g n : n £ Z} be an orthonormal basis for L 2 (0, 1) and set g n ,m = S m g n . It 
follows that {g n . m : n, m g Z} is an orthonormal basis of L 2 (R). 

Let d denote the n almost sure value of tr( / P w ) and use E to denote integration 
over a given probability space. One can then observe 



E n (tr(P w )) 



n,m / 
n,m 

= D / (Sn,^Xj(-A + V0£ mSn )di7(V) 

= E / (Sn,»(-A+%)d^). 



_ /X 

The last equality holds because v is invariant under the transformations (S r ) re m.- 
In particular, d is either or infinity according to whether or not J x (g n , Xi(~ + 
V)<7„) dP(V) is for all n or positive for some n, respectively. 

It follows that E cannot have an isolated point. To see this, observe that if E £ E 
is isolated, one can choose an interval I with /HE = {E}. But since solution spaces 
are two-dimensional, this would require 1 < tr("P w (/)) < 2 (for ii almost every ui), 
a contradiction. Hence, E lacks isolated points. □ 
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3. The Lyapunov Exponent 

We can associate a discrete SL(2, R) cocycle to our model as follows. For a given 
energy E, to £ 0, and x n ,x £ R, there is a unique unimodular matrix A E ,w(x, xq) 

so that A E ^(x,x ) ^ V y (^.°) ) = ( V y(x) ) foran y solution y 0I ^y" + K;y = Ey. 
We note that A obeys the differential equation 

^^(x^xq) = B e ^(x)A e . u] (x,xq) 



A(x ,x ) = 



1 
1 



where B E ^(x) — ^ ^ ^"(^ ^ ^ We may note that this yields an elementary 
bound on the growth of these transfer matrices: 

||4E iU ,(a;,a;o)|| < exp ( / \\B E , u (x)\\dxj = exp ( max(l, |K,(a;) - E\) Ax 

Now, for each n £ Z+, we let A E w denote the transfer matrix over the subword 

lu . . . oj n -i of lu of length n, i.e., A E cj = A E , u (£ Uo H + £ LOn _ l , 0). One readily 

sees that the cocycle condition holds: 

\n+m Am \n 

In particular, we may invoke Kingman's Subadditive Ergodic Theorem to see the 
following. 

Proposition 3.1. For each E £ R, there is a constant L^ BC (E) such that 
lim^oo i log \\A Ei J\ = L disc (E) for fj, almost every uj eft. 

This discrete cocycle's asymptotic behavior is related to the continuum cocycle's 
behavior in a natural way. 

Proposition 3.2. Enumerate A = {ai, . . . , a^} and define Q,j = {lu : luq = aj}. Let 
s = £ ai ^i(fli) + - ■ -+£ ak [i(Q,k), the weighted average of the lengths. For almost every 
lo, the continuum Lyapunov exponent L(E) = lim^oo ^ log || A EjU (x, 0)|| exists and 

satisfies L(E) = Ldl5 ; (£) . 

Proof. Consider / : £1 — > R given by f(uj) = £^ a Then, 



N ^ ' N 

n=0 



Hence, by Birkhoff 's ergodic theorem, we have 

£uq ~l~ * * * ~l~ £uj n —i 



Um = E ( /) = s 



for fj, almost every we!! 

Since lim„^ 
lu, we have that 



Since lim^oo ^ log u || and lim„^oo £ " 0+ 1 exist for (j, almost every 



,. 1 i II A n || L disc (E) 
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Thus, we see that the limit defining the continuum Lyapunov exponent exists and 
equals s~ 1 L ( ^ sc (E) along a subsequence. 

Next, we want to estimate the error in approximating the matrix A E ^(x, 0) with 
a nearby point Ab, u (j;„,0), where x n = £ Uo + ■ ■ • + £ Un _ 1 - To that end, assume 
given x > and choose n with x n < x < x n+ \. We observe the following estimates: 

- log \\A E<u (x, 0)|| - — log \\A E<u (x n , 0)|| < — log\\A E ^(x,x n )\\ 

X Xn Xfi 

and 

— log \\A E ,Mn, 0)|| - - log 0)11 

X n X 

< (— - - ) log 11^4^,0,(31, 0)11 + — log||A B ^(x„,a;)||. 

As x goes to infinity, these estimates go to zero. In particular, there is a full 
measure set of ui for which L(E) exists and is equal to Ldi SC (E)/s □ 

With this relationship established, we may consider Z = {E : L(E) = 0} = {E : 
idisc = 0}. The celebrated theorem of Ishii-Pastur-Kotani extends naturally to our 
setting: 



Proposition 3.3. S ac = Z 



-ess 



2.2 



We will 



Proof. Let fi, /i, X, D, and V^^ be defined as in the proof of Lemma 
show that v is ergodic with respect to the family of maps {S r : r G R}. 

Invariance is clear. To prove ergodicity, we follow Kirsch's argument from |29j . 
Assume that A is invariant under (>SV) re R. For each t, we may consider the section 
A t = {uj G Q : V u ,t € A}. Notice that 

T- 1 ^] V Tu ,t eA} = {coeQ: V u , t € S^ o [A] = A) = A t 

for each t, so that /J,(A t ) = or 1 for each t by ergodicity of (fl,fi,T). Moreover, 
we have 

At = {oj G n : V Utt G A} = G n : K,, G S« Wo [A] = A} = A 



and 



for all i. We then observe that 0(A) = ^(^[A]) = [ fi(A t )dt 

Jo 

hence v(A) is or 1. It follows that v is ergodic. 

We can define the Lyapunov exponent on Q by L(E,w,t) = 
lim^oo i log \\A Et ( u ^(x, 0)||, where A E ,( u ,t)(x, 0) is the natural SL(2,R) co- 
cycle associated to if Wj t = —A + V^.j. Standard arguments show that L(E,u,t) 
exists for jl almost every (uj,t) and is jl almost surely equal to a constant L(E). 
One sees that A E ^ w ^(x, 0) = A E ^(x + 1, t). 

Let Si denote the full measure set of V u ,t for which 

lim i log \\Ajs+, t t(x, 0)|| =L(E). 

x^roo X 

Next, let S2 consist of those such that V u> t G S\ for some t (this set has full 
v measure). Next, let S3 consist of those for which 

lim - log \\A E ^{x, 0)|| =L(E). 

z— >oo X 
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Since S2 and S3 have full v measure, we may choose V u 6 S% n 63 and then t such 
that V Ut t € iS'i. But then, we have 

L(E)= lim - log \\A EtU (x, 0)|| = lim - log \\A Et{Uit) (x, 0)|| = L{E). 

x— too X x— ^oo X 

By standard arguments, there is some nonrandom set £ ac such that E ac = 
v&c{Hu,t) for ft almost every (00, t). Since H u t is unitarily equivalent to H u , it 
follows that E ac = £ ac . Moreover, the Ishii-Pastur-Kotani theorem jTUJ[TT] implies 

— CSS pec O^^S 

that S ac = L- 1 (0) = Z . Thus, £ ac = Z . □ 

4. Absence of Absolutely Continuous Spectrum 

Next, we would like to prove that the absolutely continuous spectrum E ac is 
empty for aperiodic models. However, there is a slight complication - even if (Q, T) 
is aperiodic, it does not necessarily follow that the potentials are aperiodic. For 
example, we could have f a {x) = fb{x) = sin(x), £ a = 2n, £b = 47r. In this case, an 
aperiodic word containing only a's and b's will give rise to the periodic potential 
V(x) = sin(x). To remove degenerate potentials such as this, we must impose the 
following addditional assumptions: 

• Aperiodicity: The subshift and the collection {/„ : a € A} are such that 
the potentials V u are not periodic. 

• Irrcducibility: There is some I > such that the following holds. 

Suppose two "pieces" of potential satisfy X[o,£) ( fa± \\fa 2 I \fa k ) — 

X[o,£) ( /bi I/62 1 ' ' ' \fbSj (Note that this requires £ < i ai + ■ ■•4 fc ,4 I + 
• • -ib n )- Then ai = b\. In some sense, the potential has a unique decom- 
position into a concatenation of the finite pieces {/„:«£ i}. 
Notice that aperiodicity need not imply irreducibility. To see this, consider 
fa(x) = fb(x) = sin(x), f c {x) = sin(x - w), t a = 2tt, l b =l c = n. 

Under these conditions, the potentials satisfy the simple finite decomposition 
property, as described in [32] . This enables us to prove the following. 

Proposition 4.1. S ac = 

Proof. By 32 , Theorem 4. 1] , if the restriction of if u to either half line has nonempty 
absolutely continuous spectrum, then said half-line restriction must have a poten- 
tial which is eventually periodic. As the potentials under consideration are not 
eventually periodic, it follows that both half-line restrictions have empty absolutely 
continuous spectrum. By general considerations, the absolutely continuous spec- 
trum of the whole- line operator ff^ is equal to the union of the absolutely continuous 
spectra of its half line restrictions. Hence, E ac = 0, as desired. □ 

By general measure-theoretic considerations, this immediately implies 

Corollary 4.2. The Lebesgue measure of 2 is zero. 

Proof. Suppose 5 = 0. Then, given a compact interval I , since the essential 
closure of S is empty, I is covered by finitely many open intervals J\ , . . . , J m such 
that I Jfe ("I S\ = for 1 < k < m. But this implies \S fl f| = for all compact 
intervals I and hence |5| = 0. In particular, the essential closure of Z is empty, so 
Z must have zero Lebesgue measure. □ 
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5. Absence of Non-Uniform Hyperbolicity and Cantor Spectrum 

In this section we prove our main result in the general setting, namely a sufficient 
condition for zero-measure Cantor spectrum. 

Definition 5.1. Let Ae.u(x,0) be as above. We say that E E UT-L if there are 
constants C, 7 > such that \\Ae iU (x, 0)|| > Ce 7 '^' for every x £ R. That is, the 
transfer matrices grow exponentially in x and uniformly on Q,. 

Uniform hyperbolicity has several equivalent formulations; compare, for example, 
|44) . One of them is used in the proof of the following proposition. 

Proposition 5.2. S C R\UH 

Proof. Suppose E G WH. Then, for every uj £ f2, there exist linearly independent 
solutions u± of the equation —u" + V^u = Eu with L 2 decay at ±00. A standard 
argument then shows that 

_ u^(mm(x,y))u+(ma,x(x,y)) 

U_(X)U+{X) — U-(X)U + (X) 

defines the integral kernel of (H^ — E) -1 . In particular, E ^ E. □ 

Recall the following definition of condition (B) for a minimal subshift, introduced 
by Boshernitzan in [5]. 

Definition 5.3. Let ft be a minimal subshift. It satisfies the Boshernitzan condition 
(B) if there exists a T-ergodic measure \x such that 



limsup min n ■ /i ([w]) I > 0. 

n-foo \toEWfi(n) J 

Here, [w] denotes the cylinder set determined by a finite word w and yVn(n) denotes 
the set of words of length n that occur in elements offl. 

We remark that condition (B) implies unique ergodicity of (fi, T) [7]. All sub- 
shifts generated by primitive substitutions (this includes in particular the case of 
the Fibonacci substitution considered later in the paper) and, more generally, all 
linearly recurrent subshifts satisfy (B) . See [21] for many more examples of subshifts 
satisfying (B). 

Proposition 5.4. Suppose (Q,T) is a minimal subshift which satisfies (B). Then 
Z = R\UH. 

Proof. By [20, Theorem 1], it follows that ^ log \\A^, u \\ converges to L^i sc (E) 
uniformly on f2. From this and the estimates in Section 3, it follows that 
- log 0)|| converges uniformly on ft to L(E). But then for E £ K, we 

have E£Z <^ L{E) > <^ E € UK. □ 

We are now in a position to prove zero-measure Cantor spectrum. Here, we 
deviate from standard conventions slightly and call a subset of K a Cantor set if it 
is closed, does not contain any isolated points, and has empty interior. 

Corollary 5.5. Suppose (O, T) is a minimal subshift which satisfies (B) and such 
that the associated potentials V u are aperiodic and irreducible in the sense described 
in Section^ Then E is a Cantor set of Lebesgue measure zero. 



10 



D. DAMANIK, J. FILLMAN, AND A. GORODETSKI 



Proof. As it is the spectrum of an operator, E is closed. Lemma [2T2| tells us that E 
lacks isolated points. By Propositions |5.2| and |5.4| we have SCZ. The opposite 
inclusion, follows from general principles; compare, for example, the proof 



of [TTJ Theorem 2.9]. Thus, E = Z. Combining this with Corollary 4.2 which says 
that \Z\ — 0, it follows that £ is nowhere dense and of zero Lebesgue measure. 
Hence, S is a Cantor set of zero Lebesgue measure. □ 

As there are many aperiodic subshifts that satisfy (B) [3T] and most choices 
of local potential pieces yield aperiodic and irreducible potentials V u , this result 
provides a large family of continuum Schrodinger operators whose spectrum is a 
Cantor set of zero Lebesgue measure. As pointed out in the introduction, to the best 
of our knowledge, no previous examples with this spectral property were known. 



6. The Case of the Fibonacci Subshift 

In this section we study a special case, namely the subshift generated by the 
Fibonacci substitution. The alphabet is given by A — {a, b}. The Fibonacci 
substitution is the map S(a) = ab, S(b) = a. This map extends by concatenation to 
finite words over A, as well as one-sided infinite words over A. There is a unique one- 
sided infinite word u over A that obeys u = S(u). It is obtained as the limit (in the 
obvious sense) of the sequence of finite words {S n (a)} ne z, + . The Fibonacci subshift 
il C A z is defined by il — {to G A z : every subword of u is a subword of u}. It 
is easy to see (cf. [2T]) that it satisfies (B). Given a choice of i a Ab > and real- 
valued functions /„ G L 2 (0,£ a ), ft G L 2 (0,£b), we consider as above the Schrodinger 
operators {ff u } u gsi- As usual, we assume that aperiodicity and irreducibility hold 
unless otherwise stated. 

We may apply Corollary |5.5| and deduce the following: 

Theorem 6.1. There is a Cantor set SCI of zero Lebesgue measure such that 
o~(H u ) = £ for every u) G f2. 

Our goal is to go beyond this statement and to study fractal properties of S, that 
is, quantities such as the local Hausdorff dimension at a point of E. The properties 
of the spectrum are reflected by the dynamical properties of the so-called Fibonacci 
trace map. Namely, the set of all energies corresponds to a curve of initial conditions 
(which is model dependent) , and a given energy belongs to the spectrum if and only 
if the positive semiorbit of the corresponding initial condition under the action of 
the trace map is bounded. Two specific cases (piecewise constant potential and 
Kronig-Penney model) have been considered in the physics literature most often 
(mostly via numerical experiments) ; see [3J H3 1301 |3H HTJ H2] • Here we provide 
rigorous results confirming and explaining the previous numerical observations for 
each of these models. It is interesting to notice that there is an essential difference in 
the spectral properties between these models. Namely, in the Kronig-Penney model 
there are values of the energy that belong to the spectrum regardless of the value 
of coupling constant, and where the local Hausdorff dimension of the spectrum is 
equal to one (therefore the so called "pseudo bands" are formed, see also (3J). At 
the same time, these pseudo bands do not appear in the piecewise constant case, 
where the Hausdorff dimension of the spectrum in any compact domain tends to 
zero as the coupling constant tends to infinity. 
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6.1. Trace Map, Fricke-Vogt Invariant, and Local Hausdorff Dimension 
of the Spectrum. As is well known in the discrete case, the spectrum (and many 
spectral properties) of the Fibonacci model can be described in terms of the trace 
map. Let us make this connection explicit. Consider the solutions of the differential 
equation —u"(x) + f a (x)u(x) — Eu(x) for real energy E. Denote the solution 
obeying tt(0) = 0, it'(Q) = 1 (resp., it(0) = 1, u'(0) = 0) by u a>D (',E) (resp., 
w a.jv('i E)). Similarly, by replacing f a with f b , we define Ub : £,{-,E) and u b: N(-,E). 
Then, we set 

'Ua,N{l*,E) u a . D (£ a ,E) s 
.u' a , N (ta,E) u'^ D (e a ,E) / 



M(a, E) 

M(h w\ - f u b,N(h,E) u bjD (e b ,E)\ 
M ^^-\u' biN (£ b ,E) u' bjD (e b ,E)J 



and 

(1) x- 1 (E) = ±tr(M(b,E)), 

(2) x (E) = ^tr(M(a,E)), 

(3) Xx(E) = itr {M(b, E)M(a, E)) . 

The map E h-> (x%(E), X\(E), Xq(E)) will be called the curve of initial conditions. 
The trace map T : R 3 — > M 3 is given by 

T(x,y,z) = (2xy - z,x,y). 

Define 

I(x, y, z) = x 2 + y 2 + z 2 — 2xyz — 1. 

Again, as is well known, the trace map T preserves / and hence its level surfaces 

Sj = {(x,y,z) EM. 3 : I(x,y,z) = I}. 

In particular, all the points T n (x\(E), Xq(E), x_i(E)) lie on the surface Sj^e)i 
where (with some abuse of notation) we set 

I{E)=I{x 1 {E),x {E),x- 1 {E)). 

The surfaces Sj undergo a transition at / = 0. For negative values of /, Si has 
one bounded connected component and four unbounded connected components. 
For 1 = 0, the four "outside" components attach to the inner part in four conic 
singularities. These singularities resolve as I becomes positive, and in this case 5/ 
is connected and smooth. We show plots of Si for I < 0, I = 0, and / > in 
Figure 1. 



We will need the following statement about the dynamics of T on Si for / < 0. 

Lemma 6.2. Suppose I < 0. Then, every point on the bounded connected com- 
ponent of Si has a bounded T-orbit, and every point on one of the unbounded 
connected components of Si has unbounded forward and backward orbit under T. 

Proof. By ^ Proposition 4], the bounded connected component of Si is forward 
and backward invariant under T and hence the T-orbit of any point on it is bounded. 

On the other hand, by [381 Theorem 4.3], any point on one of the unbounded 
components of Si has unbounded forward and backward orbits under T. □ 
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Figure 1. The surface Si for / < 0, I = 0, and I > 0. 



The dynamical spectrum B is defined by 

B = {E e R : {T"(zi(£),x (.E),a:_i(.E))} neZ+ is bounded}. 
Proposition 6.3. We /iave S = B. 

Proof. Recall that the one-sided infinite Fibonacci word u is the limit of the finite 
words S n (a), n > 0. Moreover, it is not hard to see that the Fibonacci subshift 
Q contains an element oj u whose restriction to the right half-line coincides with u. 
Since the spectrum of H u is w-independent by minimality, we can study the set 
£ by considering the spectrum of the particular operator H Uu . For E £ R and 
n > 0, denote by M n (E) the transfer matrix corresponding to the operator H Uu , 
the energy E, and the interval arising from the word S n (a) through our substitution 
procedure. Due to S n+1 (a) = S n (ab) = S n (a)S n - 1 (a), we have 

(4) M n+1 (E) = M n -i(E)M n (E), n > 1. 

Using the Cayley-Hamilton theorem, one sees that for x n (E) :— ^Tr M n (E), n > 0, 
we have 

(5) x n+1 {E) = 2x n {E)x n - X {E) - x n - 2 {E), n>2. 

Note that Q and ^ can be used to define M n (E) and x n (E) for n < 0. In 
particular, 

M- X {E) = M 1 (E)M (E)- 1 . 

Since M\{E) corresponds to S 1 (a) = ab and Mq(E) corresponds to S°(a) — a, this 
shows that the matrix M_i(E) is just the transfer matrix across ff, corresponding 
to the energy E. In other words, our current definition of x n (E) for n — —1,0, 1 
agrees with the definition given above; compare Q-Q. 
In particular, we have 

B = {E e R : {x n {E)} n&+ is bounded}. 

We can now prove the inclusion U B C £." Suppose E € B. As is well known 
(and in fact easy to check), u does not only have S n (a), n > as prefixes, but 
also S n (a)S n (a), n > 2. Combining this with the boundedness of x n (E) for E e 
B and the Gordon lemma (see, e.g., [TH |2H [2SJ)> we find that no solution of 
— u" + V UJu u = Eu is square-integrable at +oo. This implies by standard arguments 
that E £a(H Uu ) = S. 

Next, we prove the inclusion "S C B." Denote 

<j n = {E eR:\x n {E)\<l}. 
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By Floquet-Bloch theory, a n is the spectrum of the Schrodinger operator H n with 
periodic potential obtained by repeating the piece corresponding to S n (a). The 
operators H n converge in strong resolvent sense to H u for some suitably chosen 
u G Q. This implies 



(6) E = a(H u ) C P| |J a k . 

n>0 k>n 

On the other hand, a minor modification of (37J Proposition 12.8.6] shows that 
E G" B if and only if there exists n with |x„+i(£')| > 1, |x n (J5)| > 1, and 
\x n+ i(E)x n {E)\ > \x n -i(E)\. Moreover, in this case, we necessarily have |xfe(-E)| > 
1 for all k > n. 

Combining these two facts, we can now conclude the proof. Suppose E B. 
Then there exists n with \x n+1 (E)\ > 1, |a;„(i?)| > 1, and \x n+1 (E)x n (E)\ > 
\x n -i(E)\. Choose e > such that \x n+ i(E')\ > 1, \x n (E')\ > 1, and 
\x n+1 (E')x n (E')\ > \x n -x{E')\ for every E' with \E - E'\ < e. For each of these 
E' , we therefore have \xk(E')\ > 1 for all k > n. This shows that 



k > n 

and hence E g S by (|6f. □ 

At this point we would like to point out a major difference between the discrete 
case and the continuum case. In the discrete case, 1(E) is actually constant, while 
in the continuum case, it is in general not constant (as we will explicitly see below). 
As we will see later, this leads to new phenomena in the continuum case that make 
its study worthwhile. 

Moreover, in the discrete case, the invariant is always non-negative. In the 
continuum case, we cannot a priori rule out that it may be negative for some 
energies. The following proposition shows that even if this happens, these energies 
are not essential as they must lie outside the spectrum. As a consequence, the 
study of the dynamics of the trace map on Si for I > 0, which has been investigated 
heavily in the papers on the discrete Fibonacci Hamiltonian, is sufficient to describe 
the corresponding spectral properties of the continuum Fibonacci Hamiltonian. 

Proposition 6.4. We have 1(E) > for every fieE. 

Proof. Assume there is E £ E with 1(E) < 0. Then, since the forward trace map or- 



bit of (x2(E), xi(E), xq(E)) remains bounded by Proposition 6.3 this point must 
belong to the compact component of the invariant surface Sj^e) by Lemma 6.2 
Now wiggle E. By continuity of x n (-), nearby E"s must have I(E') < anc 
(x2(E'),Xi(E'),Xo(E')) contained in the bounded component of Sj^e) as well. 



Thus, again by Lemma 6.2 (xi(E'), X2(E'), x^(E')) has bounded forward trace 
map orbit, too, and hence E' S E by Proposition |6.3| It follows that E is not a 
Cantor set, which contradicts Theorem |6.1| □ 



We can now address the local fractal dimension of the spectrum. 

Theorem 6.5. There is a continuous map D : [0, oo) — > (0, 1] with the following 
properties: 

(i) dim l £ c (E, E) = D(I(E)) for every E G E. 

(ii) We have D(0) = 1 and 1 - D(I) ^Vl as 1 10. 
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(iii) We have 



lim D(I) •log/ = 21og(l 

I— >oo 



(iv) D is real analytic in (0, oo). 
Proof. Reca ll th at by Lemma |2.2[ £ contains no isolated points. Therefore, by 



Proposition [^i] and [23] Theorem 2.13], dim^ c (E,j:) for EeS depends only on 
the value of 1(E). This implies (i). Given this, (ii) follows from [17] and (iii) 
follows from [T5]. (iv) follows from Theorem 5.23]. Together with (ii), this 
implies continuity of D on [0, 00). □ 

6.2. Explicit Computations and Formulae in Special Cases. In this sub- 
section we consider special choices of the local potential pieces f a , fb for which it 
is possible to compute the Fricke-Vogt invariant 1(E) explicitly. As we have seen 
above, this has a direct connection to the local structure of the spectrum near an 
energy The explicit formulae also permit us to identify the limit behavior 

of the local Hausdorff dimension of the spectrum in the regime of large energy or 
small/large coupling. 



6.2.1. A Piecewise Constant Potential. For comparison purposes, let us start with 

. d 2 

dx 2 



the free case, f a = fb = 0. That is, we consider the free Schrodinger operator 



in L 2 (M.), but view its (zero) potential as tiled according to a Fibonacci sequence. 
In this degenerate case we of course do not have aperiodicity and irreducibility. 
However, all quantities and formulae that arise from the presence of Fibonacci 
symmetries still exist. In particular, we may compute the curve of initial conditions 
and the resulting Fricke-Vogt invariant 1(E). 
For E > 0, we have 



/ cos \J~E~ — L sin \J~E~ 
M(a,E)=M(b,E)=[ y v 

\-VEsmVE cosVE 



and hence x-i(E) = xq(E) — cosvE. Similarly we obtain x\(E) = cos(2\/E). 
For E < 0, we have 

/ cosh V— E A-r sinh \l — E\ 

M(a, E) = M(b, E) = , , J_ . 

V ' ' V ' ' ^y^Bsinhy^B coshV-B J ' 

and hence X-i(E) = x (E) = cosh \J—E, x\(E) = cosh(2V— E). 

A direct calculation shows now that in both cases 1(E) = 0. By continuity, 
we have also 1(0) = 0. Thus, as in the discrete setting, the invariant vanishes 
identically in the free case. This suggests that the dynamical behavior of the trace 
map near So is essential for a study of the continuum Fibonacci Hamiltonian in the 
small coupling regime (at least in compact energy regions), just as it was the case in 
the discrete setting and which has led to numerous recent advances [TBI H3 HH1 US] • 



Consider now the case f a = A-X[o,i) and fb = 0-X[o,i)> where A > 0. Clearly, when 
A = 0, we obtain the free case considered above; and when A > 0, the resulting 
potentials are aperiodic and irreducible. Let us compute the initial conditions 
(xi(E), x (E), X-i(E)) first in the case E > A. We have (compare with (2.15) from 
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my 

x_i(E) = cos VE, 
x (E) = cos VE — A, 

x 1 (E) — cos VE cos \/E — A — ~ {^J X \/ ^E ^ S ™ S ™ ~ ^' 
From this we can calculate /(E) explicitly: 
(7) 1(25) = 1 — ^- - sin 2 VE sin 2 Ve^A. 



Notice that if sinVE = 0, then cos VE e {-1,1}, and we have either 
(xi(E),xq(E),x-i(E)) — (cos \JE - A, cos \JE - A, 1), which is close to (1,1,1) 
if -E 1 > 1, or (xi(E),x Q (E),x^i(E)) = (- cos VE - A, cos VE - A, -1), which is 
close to (1,-1,-1) if E > 1. 

Let us now calculate the initial conditions and 1(E) for E g (0,A). We have 
(compare with (2.14) from [41]): 

x_ x (E) = cos vE, 
x (E) = cosh VE — A, 



xi (E) — cos VE cosh VA — E + i f y - — y sin VE sinh VA - E, 

and hence (compare with (3.4) from |41j ) 
(8) /(E) = ^ g(A A ' g) sin 2 VEsinh 2 VA~E. 

With these explicit formulae, we can now establish the asymptotics of the local 
Hausdorff dimension of the spectrum in the small coupling regime as well as in the 
high energy regime. 

Corollary 6.6. We have the following asymptotics: 

lim inf dimj? c (E, S) = 1, 
lim inf diml? c (E,S) = 1. 

K^oo £6£n[if,oo) 

Proof. From Q and ^ we know that 1(E) -> as A -> 0, uniformly in E > 
(hence uniformly in E G £). This implies the first statement of Corollary 6.6 due 
to Theorem 16.51 

Also, from ^ we see that for any fixed value of A, we have 1(E) — > as E — > oo, 
which implies the second statement, again by Theorem |6.5| □ 



In the large coupling regime, we have the following result: 
Corollary 6.7. For any compact 5cR, we have 

lim dim H (Sn S) = 0. 

A— yoo 
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Proof. Consider E e S with 1(E) < N 2 , and so that A > E,N. Then, due to 



\Jl(E) ~ A 1/2 sin VE sinh s/X- E 



<N, 



and hence |sin\/l£| < "T/a^T 1 ^ ^" This implies that jcosV-Ej > 1/2. Now we 



have 

and 
\ Xl (E)\ 



!(E)\ = Icosv^l < 1, \x Q (E)\ = |coshVA-£;| - e Al/2 > 1 



EcoshVX- E 



1 



1 X — E 



> 



,1/2 



> e AV2 

> 1. 



2 \ V E 

CX 1/2 sin Vi? sinh VX - E 
N 



E 
E-X 



sin \AE sinh VX-E 



This implies £ ^ B by [351 Lemma 2] and hence ^ E by Proposition 6.3 
Therefore, for any E £ S H E, we must have > N 2 if A is sufficiently large, 

and this proves Corollary |6.7| □ 



One in general expects that in the regime of small coupling or high energy, the 
local characteristics of the spectrum should approach those of the free case. Simi- 
larly, one expects that for fixed energy (or fixed compact energy region) , the poten- 
tial should dominate in the large coupling regime and hence the situation should 
become as singular as possible. Corollaries |6.6| and |6.7| show that the expected be- 
havior holds true in the simple situation of a piecewise constant potential, where the 
Fricke-Vogt invariant can be computed explicitly. It is quite natural to expect that 
these results should extend to much more general choices of the potential pieces. 
In fact, we ask the following: 



Question 6.8. Is it true that Corollaries \6.6\ and \6. 7| hold regardless of the shape 
of the bump? That is, if we replace f a = A ■ Xfo. i) a n d fb = ■ X[o,i) by general 



f a € L (0,£ a ) and f e L (0,£ o ), do Corollaries 6.6 and 6.1 continue to hold as 
stated? 



6.2.2. The Kronig-Penney Model. The Kronig-Penney model places local point in- 
teractions at a discrete set of points. Here we are interested in the case where the 
location of these points is dictated by a Fibonacci sequence. This model was con- 
sidered in [5] in an "off-diagonal" setting. Let us make the calculations also in the 
following ("diagonal") setting: We take £ a — lb — 1 and f a (x) = XS(x), fb(%) = 0. 
Note that strictly speaking, this model is not contained in the general class of mod- 
els considered earlier in the paper. We nevertheless compute the invariant in this 
case as, together with our earlier results, this is quite constructive. We leave the 
extension of the general part to the interested reader (note that [32] actually did 
consider measures rather than potentials). 
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For E > 0, we have 
M(b,E) 




( cos \J~E~ -i= sin \/ E 



cos + J= sin V^E 
A cos \J~E~ - V~E sin \[E 




In this case we have: 

x-i(E) — cosy/E, 
x (E) = cos 2y/E + sin y/E, 

xi (E) = cos 2-/E + ^= sin 2-/B. 
Explicit calculations show that in this case 

(9) 1(E) = ^ sin 2 VE 

(compare with Section 2 in [3]). Notice that, irrespectively of the value of the 
coupling constant A, if sin^/E = for some energy, then the formulae above show 
that this energy is in the spectrum, and the corresponding point is on internal 
part of the Cayley cubic. Since there are no isolated points in the spectrum, this 
implies that at those points the local Hausdorff dimension of the spectrum is equal 
to one (which explains the nature of "pseudo bands" mentioned in [3] , see also [M] 
where a similar observation was made) , and in all other points of the spectrum, its 
local Hausdorff dimension is strictly less than one. In particular this shows that 
the analog of Corollary |6.7| does not hold in this case, which shows an essential 
difference between the two simple models we have considered. 

Certainly, due to ([9]), the analog of Corollary 6.6 does hold for this model. This 



shows that at least the heuristics leading to the general expectation in the small 
coupling and high energy regimes do apply in this case. 
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